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1. Introduction 

In ordinary zero-temperature perturbative field theory, this conference series has witnessed a 
rapid development of automated computer-algebra aided approaches over the past decade or so. 
Triggered perhaps by the integration-by-parts (IBP) method ^ ^ that received an enormous boost 
from its algorithmic description by Laporta [^], the problem of reducing the large variety of Feyn- 
man integrals that occur in a given perturbative computation to a small number of master integrals 
can nowadays be regarded as solved, in the sense that there exist independent public computer- 
algebraic implementations [Q-|6|] which, given sufficient hardware resources, perform this step in 
an automated way. Concerning the subsequent evaluation of master integrals, the degree of au- 
tomation is less developed - after all, integration is still an art. However, besides a much better 
understanding of the analytic structure of the numbers and functions that are contained in the per- 
turbative series [|7|-[I(J, there has been progress in developing powerful approaches to obtain ana- 



lytic results ^ [IT]] as well as stable and fast numerical tools [12, 13] that promise to be applicable 
to large classes of Feynman integrals. 

In this note, we will discuss the corresponding situation for finite-temperature perturbation 
theory. The present state of affairs is that, while many of the methods and algorithmic tools devel- 
oped for zero-temperature field theory (such as generation and classification of Feynman graphs, 
efficient color and Lorentz algebra as well as IBP reduction methods) can be - and have been 
- applied with only minor adjustments to finite-temperature systems, only very few master sum- 
integrals that remain after the IBP reduction step are known (see, e.g. Refs. |[T^-|r^] or the review 



[]19[]). It seems difficult to profit from the comparably mature zero-temperature techniques for this 
step, due to the very different analytical structure that the sums bring about. Here, we will there- 
fore give only a brief discussion on thermal IBP methods, to concentrate then on sum-integral 
evaluation, producing a new result for a phenomenologically relevant case. 

2. Sum-integral reduction via IBP 

The basic principle of the IBP method [|l|, applies to any integral, and can hence also be 
applied to the sum-integrals that occur in finite-temperature field theories. The Matsubara sums 
are simply left untouched, while their summands are interpreted as massive loop integrals in a 
reduced space-time dimension and with masses provided by the Matsubara frequencies, with the 
IBP relations acting upon these massive loop integrals. A more detailed review of these techniques 



is given e.g. in Ref. Q20|]. Similar to the zero-temperature case, the linear relations among sum- 
integrals - as generated by IBP on the Matsubara summands - can be systematically used for a 
sum-integral reduction step, employing a variant of the Laporta algorithm [g]]. This approach has 
been successfully used for a number of higher-order calculations in finite-temperature QCD (for 



examples, see e.g. [EQ- 22]), to enable basis transforms among master sum-integrals (see below), 



or even to aid in managing the infrared (IR) behavior of Matsubara summands when evaluating 



master integrals [23]. To show two concrete examples, IBP has revealed that a non-trivial 2-loop 



sunset-sum-integral [ |24| ] and a non-trivial 3-loop mercedes-sum-integral vanish identically: 



i 

Tp 



pqP 2 Q 2 (p-Q) 2 



^0, jTnn =0, (2.1) 
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where we have used a generic notation for massless bosonic 3-loop vacuum sum-integrals 

jafiy J (Pb) a (Qo) P (Ro) r (22 s 

abcdef h Q R[P 2 nQ 2 } b [R 2 Y[(P-Q) 2 } d [(P-R) 2 Y[(Q-R) 2 y ' ^ ' ; 

In our notation, bosonic (Euclidean) four-momenta are denoted by P = (Po,p) = (2nnT,p), with 
T being temperature of the thermal system, inverse propagators are P 2 = P^ + p 2 , and the sum- 
integral symbol stands for 



JP ne7L J 



d rf p 

with d = 3-2e . (2.3) 



{2n) d ' 



A class of 1-loop bosonic tadpoles that we will need below can be evaluated analytically as 

r i _ 2 T q2 S -d) r(5-f) 

s ~Tq[Q 2 Y {2KTf>- d (An) d / 2 T(s)' y '' 



In Ref. [22], it has been shown that the computation of NNLO corrections to the spatial string 
tension of pure Yang-Mills theory (mapped to Taylor coefficients of background-gauge-field self- 
energies [21]) can be reduced to 3-loop basketball-type master sum-integrals Kj (and products of 
simpler 1-loop ones /,•), giving the gauge-invariant expression 

/ 9 \ 



HtsW =C\ IVPiWKi + PioWlM+PnWhhh+PnWhhhj (2.5) 

j F \ / j-000 tOOO t002 t020 t220 t400 t800 ,730 r820 \ ,n ^ 

•[Ai, . . . ,A9j = V22001l5- / 31001l5- / 320011i- / 410011i- / 510011'' / 510011i- / 710()lli- / 730-llli- / 82O-lll/ ) \ z - ) 



where the coefficients j8, are rational functions that contain simple and double poles in l/(d — 3), 
and where all massless bosonic 3-loop vacuum sum-integrals Kj have been represented in terms 
of the notation introduced in Eq. (|2T2|). After an extensive reverse search in our IBP database, it 
proves possible to transform the above expression into the equivalent representation 

Ht 3 (0) = C 3 A ^£ n(d) Vi + r 7 (d)hhh + n{d)hh_h + r 9 {d)hhl^j (2.7) 

r T/ x/ l _ r tOOO 7000 7002 7020 7022 ,000 1 n o\ 

\yU- ■ ■ ) v 6f = \ , '121110' 21 11 10 '"'221 110 ' 31 11 10 ' 41 11 10 ' 310011 J 



which, much in the spirit of the epsilon-finite basis advocated in Ref. []25|], does not contain diver- 
gences in the coefficients n(d) as d — > 3 (note that Ve = K2 was already contained in the old basis 
listed in Eq. (2.6), and was kept since it is has already been evaluated in Ref. JT8|]). Eq. (2.7) can 
now serve as a convenient starting point for determining the spatial string tension, once the five 
unknown master sum-integrals {V\,. .. ,Vs} have been evaluated up to their constant parts. For V\, 
see the next section. 



3. A 3-loop master sum-integral of dimension zero 

Let us now turn to a concrete example of 3-loop sum-integrals, to demonstrate the main tech- 
niques that are used in the master integral evaluation step. Here, we work in d = 3 — 2s spatial 
dimensions, and wish to evaluate the e-expansion of V\ given in Eq. ( |2T8| ) up to the constant term. 
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Vi is a sum-integral of spectacles-type, for which a generic evaluation procedure has been 
discussed in Ref. Pq]. Subtracting subdivergences, the sum-integral is split into its finite part 
(which can typically be computed only numerically in configuration space), and its divergent part 
(which can be expressed analytically in terms of Zeta and Gamma functions), treating the Po = 
mode separately (softening the IR via IBP relations if needed). We thus decompose V\ as 

Vi =J^m =^U 21 (P)U n (P) (3.1) 

= ^l{(n 21 -nf 1 )(n 11 -nf 1 )+nf 1 (n 11 -nf 1 ) + (n 21 -nf 1 -n^ 1 )(nf 1 -nf 1 )} + 
+^{(n 2 i -nM +nf 1 (nf 1 -nfr+n&n?,} +^3^ , 0.2) 

where the first line collects the finite pieces (the primed sum excludes the Po = term), the second 
line the divergent pieces as well as the zero-mode, and we have used the 1-loop 2-point structures 

_J 1 b _ [ 1 G(a,b,d + 1) 

ab( P ) = ab=J Q [Q2]a[(p-Q)2]b ' U " b = J Q [Q2]a [(/> _ Q)2]i = [p2] fl +fc-(d+l)/2 ' (33) 

where G(ji,j 2 ,<i) = /,l w/2r / Tp? uv3 T ' (3,4) 

(4n) d ' 2 r(si)r(s2)r(d - s { -s 2 ) 



L h _n [P 2 1 £ 



n c = rr B i i n D 



<* - a* [p 2 ]fc [p2] fl ' a* - ( 27 i:r) 2e ' " afc ~ t e [e 2 ] u [(p - e) 2 f 



nf fe , n* =£s e 1 (3.5) 



We will in the following evaluate the various contributions to V\ of Eq. ( |3.2[ ) in turn. 



3.1 Finite pieces 



Using inverse 3d Fourier transforms of the 2-point functions, the first term of Eq. (3.2) reads 



^ x ( n 2i- n fi)x( n n- n n) 

^7(27T) 3 P2 2(47T)2i f f Q \Q Q \ 4 J S 2 \ S J 

2 r°° roo i / i \ _-(|n|+|m|+|m+n|)v , * 

/* °° 1 / 1 \ r 

= WfL^Ty { COthy -y) [^<l-e- 2y )+y] - ^[ln(l - e^) +4y] - Uy* + 

+ 61n(l-e- 2v )Li 2 (e- 2 >') + 12j[Li 2 (\/{\-e- 2y ))-iK\n(\-e- Zy )\ -6Li 3 (l -e 2 ?) 



(4k) 



+ &(£) 

+ 0(e), ci « 0.6864720593640618954(1) , (3.6) 



where we have used dimensionless variables such as r = 27iT\r\ = x and <2o = <2o/(2ttP) = m; 
performed the momentum integration as well as the angular integration in configuration space in 
the third line; integrated over x and performed both sums in the fourth line; and performed the last 
step simply via Mathematica's numerical integration routine. 
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In full analogy 1 , we get for the second term of Eq. (3.2) 
f, 1 



^> xn fi x ( n n- n n) 



d 3 p 1 



I. Ah 



n 



-- )e- 2nx + @(e) 



1 



{An) 



6 2 3 



+ 0(e) 



(3.7) 



where G « 1.28243 is the Glaisher constant. Similarly, the third term of Eq. < J3-2| > evaluates as 
f, 1 



with / 2 io(jc,n) = e 2 " v B(^,«+l,0) +#„ - ln(l -<? _2 *) + e 2,lx Ei(-2ra) +ln — - 7e - 

n on 



\aP 2 

""87T 2 

— 2x\ i _2nr 



-2nx 



Jo 

C2 



e 2n *Ei(-2ra;)+ln^ + <)fe /2io(jc,«) + ^(e) 



2?2 



(47l) f 



+ ^(e), c 2 « -6.404(1) , 



(3.8) 
(3.9) 



where we have used the incomplete Beta function B(z,n + 1,0) = ffidtt" /(I —t) = — ln(l — z) — 
Y^n=\Z m /m, harmonic numbers H n = Y? m =\ ' l m an d exponential integral Ei(z) = J^dte 1 jt. After 
integrating over p the angular r integration was trivial, leaving Eq. (3J5) where the summation 
converges somewhat slowly and the evaluation of the integrand itself is costly since it contains 
special functions; for the numerical precision given above, we have truncated the sum at « max = 
7000 and estimated the remainder by fitting to a power-law a/n b in the interval n € [7000, 19000] 
(obtaining b ss 1 .9) and summing this fit to infinity. 

3.2 Divergent pieces 



Introducing the following 2-loop vacuum sum-integrals (see e.g. Appendix A and B of Ref. ||26|]) 



A(S[,S 2 ,S 3 



\S4\d) = % 
Tp 



s Qo \p \ S4 



L(si,s 2 ,S3;s4,s 5 ;d) 
L(211;00;d) 



pq [Q 2 Mp 2 } S2 [(p-Q) 2 Y 3 
2T 2 C(2 Sl23 -2d- s 4 ) rjsn - f )r> 12 - f )r(f - si)r(si23 - d) 

(2KT) 2s ^- 2d - s * (4K) d r{s 2 )T(s 3 )r{d/2)r(sn23-d) 

(PoY 4 (QoY 5 



(3.10) 



pq[p 2 }^[Q 2 Y 2 [(p+Q) 2 } S3 

IBP 1 



(d-S){d-2) 



hh 



(3.11) 
(3.12) 



'Here, we use (d — 3) 2P ° pi ] ^ Tlf l — jp^rH^ — 4 ^Jp°|i ■ where we have treated \Pq\ as a mass in the 3d 

1-loop self-energy integrals IPj, , and noticed that IT^ is finite at d = 3, and that the tadpole Uq 2 is known analytically. 
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and noting from Eqs. (fO|),(|3~3|) that n D ~ Yl B ~ l/[P 2 ] x are simple powers, the fourth, fifth and 
sixth terms of Eq. (32) are effectively 2-loop structures and can be evaluated analytically as 

r, i 



^(n 21 n E 2l )n n (27rr)3 _ f/ 

^ n& (nf ! - nfo = g(i, i, d + 1) [g(2, 1,^ + 1 )/ 5 - d + / 2 /( 7 - rf )/2 + /i/( 9 

jf'^n^nfj =G(l,l,rf+l)A(2,(5-rf)/2,l;0;rf) , 



L(211;00;rf) -A(121;0;d) -A(211;0;<Q 

d)/2 



(3.13) 
(3.14) 
(3.15) 



where we have used the abbreviation / v = 7 V — (2%T) d 3 I s+ (d-3)/2- All three expressions contain 
poles as d — > 3; the first two contribute starting from 1/e 3 , the last one merely from l/e. 

3.3 Zero mode 



The seventh term of Eq. ( p.2| ) can be decomposed into finite (first line) and divergent parts as 

jf tpoji n 2i n n = ^ Po ^ (n 21 - nfo (n n - n* - nf t ) + 

+^5p ^{nf 1 n 11 +n 21 (nf 1 +nf 1 )-nf 1 (nf 1 +nf 1 )} . (3.16) 

In full analogy to Sec. Oj, the first (finite) term is treated in 3d coordinate space, as 

^s Po l x (n 21 -nfj) x (n u -nfi -nfo 

4 



OO /> CO 



(47r) 6 Jo Jo 
1 4 



dr/dy X+y " |x " yl Ml- e -^ 



(47T) 6 3 

C3 



jo y V 



cothj 1 



cothj 



1 +^(£) 



4y 3 - 2% 2 y + 3 [Li 3 (e 2y ) + 27n> 2 ] - 3£(3) 



+ 



(An) 



+ 0(e) , c 3 « 10.33244698246374834(1) 



(3.17) 



In full analogy to Sec. p.2| , the second (divergent) part of Eq. ( p.l6[ ) contains 1- and 2-loop structures 
G and A only and is hence known fully analytically in terms of Gamma and Zeta functions: 

^Spoji { n fi n n + n 2infi +n 21 nf 1 -nf, (nf, +nf 1 )} 

= TG(2,l,d)A(4-d/2,l,l;0;d) + TG(l,\,d)A(3-d/2,2,l;0;d) + 

+ G(l,l,J+l)A((5-J)/2,2,l;0;J)-0 scaMr ee • (3.18) 

3.4 Result 



Collecting from Eqs. (p.6D,(P.7D,(p.9D,(P.13D, d3.14|),(|3.15D,(|3.17D,([3.18D and expanding around 
d = 3 — 2e, we finally obtain 



1 6(4?r) 6 \4kT 2 ) 



3e 



3 1 

e J e 1 e 



^ + ^ + ^(l3-6^ 2 + ^-12 7l -3C(3)) + 



37T 2 

T 
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+ (51- 42ye 2 + An 2 + ln(27l) - 121n g) + 21n(2) ( 12 - 12^ 2 - 24 7l - £ (3)) + 
+ 6^ (3C(3) - 4 - 471 ) - 84 7l - 36 72 + ^ C(3) - 16£'(3) + 6(a + c 2 + c 3 ) ) + 0(e) 



1 / 1 \ 3e / 1 2.86143 15.2646 



+ — +28.56(1) + 0(e) ), (3.19) 



6(4tt) 6 \T 2 J \e 3 e 2 

with (ci +C2 + C3) « 4.615(1), where the numerical error is dominated by Eq. (^^). The analytic 
part of the result contains the Glaisher constant G, for which 121n(G) = 1 + C,' (—1) / C,(—\), zeta 
values as well as the Stieltjes constants 7, defined by £(l + e) = l/e + ^- 7 i£ + 7 2£ 2 /2 + ^ > (e 3 ). 
Note that our 3-loop sum-integral V\ does contain a 1/e 3 divergence, as would be naturally 



expected from the fact that the 1-loop tadpoles of Eq. ( |2.4[ ) diverge as 1/e at most (recall that 
d = 3 — 2e in our notation). This is in fact the first example where we observe such behavior - all 
previously known non-trivial 3-loop cases in fact started at order 1/e 2 . 



4. Conclusions 

While automated methods for Feynman integral reduction work well when applied to the sum- 
integrals that arise in finite temperature perturbation theory, the step of evaluating the resulting 
master sum-integrals is in a much less mature state. Only a small number of non-trivial higher-loop 
sum-integrals are known so far, their evaluation resting on a case-by-case analysis with intricate 
subdivergence subtraction techniques, such as demonstrated here on the example of a new bosonic 
3-loop tadpole that enters in a NNLO determination of the spatial string tension. Typical results 
are partly analytic and partly numeric, since often the constant terms of the Laurent series cannot 
be obtained in closed form, but mapped onto simple finite low-dimensional representations. 

Future progress in the field of classification and evaluation of sum-integrals would certainly 
give a boost to the field of finite temperature field theory. It will be interesting to see whether 
this progress originates again from a fruitful application of zero-temperature methods (such as 
e.g. Mellin-Barnes representations, classes of multiple nested sums and integrals such as harmonic 
polylogarithms, or systematic numerical methods based on sector decomposition), or whether com- 
pletely new ideas and structures are needed. 

The work of I.G. is supported by the Deutsche Forschungsgemeinschaft (DFG) under grant 
no. GRK 881. Y.S. is supported by the Heisenberg program DFG, contract no. SCHR 993/1. 
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